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Abstract. The problem of factorization and parametrization of compactly 
supported biorthogonal wavelets was reduced to that of pseudoidentity ma- 
trix pairs by ResnikofF, Tian, and Wells in their 2001 paper. Based on 
a conjecture on the pseudoidentity matrix pairs of rank 2 stated in the 
same paper, they proved a theorem which gives a complete factorization 
result for rank 2 compactly supported biorthogonal wavelets. In this pa- 
per, we first provide examples to show that the conjecture is not true, then 
we prove a factorization theorem for pseudoidentity matrix pairs of rank 
TO > 2. As a consequence, our result shows that a slightly modified version 
of the factorization theorem in the rank 2 case given by Resnikoff, Tian, 
and Wells holds. We also provide a concrete constructive method for the 
rank 2 case which is determined by applying the Euclidean algorithm to 
two polynomials. 
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torization, Vanishing Moment 
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1. Introduction 

It is well-known that compactly supported orthogonal wavelets cannot be 
symmetric except for Haar wavelets O [Til [13]. In signal processing, the sym- 
metry property corresponds to the linear phase condition, and since biorthogo- 
nal wavelets permit symmetry, they are preferred over the orthogonal wavelets 
in applications such as digital filter banks. Because of the flexibility of the 
biorthogonal condition [21 |3l HJ [5l [6l [TJ |9l [13] , one can add other conditions to 
the definition of a biorthogonal wavelet according to the needs. To guarantee 
the existence of scaling functions and wavelet functions, a zeroth-order van- 
ishing moment condition was imposed on the biorthogonal wavelet matrix in 
addition to the usual quadratic condition in [11]. This extra requirement must 
be taken into account when analyzing the structure of these wavelets. In par- 
ticular, the factorization results from the earlier investigations [2l[3l[6ll9l[T0l[T2] 
do not apply to these wavelets. It was proved in [llj that any biorthogonal 
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wavelet matrix pair can be decomposed into four components: an orthogonal 
component, a pseudoidentity matrix pair (defined in [H]), an invertible ma- 
trix, and a constant matrix. Since the factorization and parametrization of the 
orthogonal wavelet matrices is well understood [HI [121 [13], the factorization 
of a biorthogonal wavelet matrix pair reduces to that of a pseudoidentity ma- 
trix pair. To obtain a complete factorization and therefore a parametrization 
of the biorthogonal wavelets, a conjecture was proposed for the pseudoiden- 
tity matrix pairs of rank 2 in ^T\, and based on this conjecture, a complete 
factorization theorem was proved for the rank 2 case therein. 

In this paper, we prove a factorization theorem for a pseudoidentity matrix 
pair of rank m > 2, which in particular implies that the factorization theorem 
for the rank 2 case stated in [11] is true (with a slightly modified statement) 
in spite of the invalidness of the conjecture proposed there. We first provide 
examples to show that the conjecture is not true and then proceed to prove a 
factorization theorem. Our proof is elementary and constructive: it applies to 
any rank and provides an algorithm for the decomposition. 

We organize our presentation as follows. In section 2, we first recall the 
definition and the basic properties of biorthogonal wavelets as given in [TT] . 
then we discuss some counterexamples to the conjecture stated in [TT]. In 
Section 3, we prove our factorization result on pseudoidentity matrix pairs and 
thus obtain a complete factorization for the biorthogonal wavelet matrix pairs. 
In section 4, we consider the rank 2 case in some detail and show that in this 
case, the factorization of a pseudoidentity matrix pair is basically determined 
by the Euclidean algorithm for finding the greatest common divisor (GCD) of 
two polynomials. We conclude the paper in section 5 with some discussion. 



2. Biorthogonal wavelets 

Let C be the set of complex numbers, let F C C be a subfield invariant 
under the complex conjugation, and let R = ¥[z, z~'^] be the ring of Laurent 
polynomials over F. The reason for us to work with F instead of just C is 
that it will allow more flexibility in the applications, since the discussions here 
apply to a whole range of fields including the rational numbers Q, the real 
numbers M, the complex numbers C, as well as other algebraic extensions of 
Q. For instance, in some applications of digital filters, it is desirable to work 
with rational numbers, since it is easier to implement the filters. We denote 
by gl{m, R) the set of all m x m matrices over R and denote by GL{m, R) the 
set of all invertible m x m matrices over R. We also denote by gl{m,¥) and 
GL{m, F) the set of m x m matrices and the set of invertible m x m matrices 
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over F, respectively. For a matrix A{z) G gl{m, R), we can write 

(2.1) A(z) = ^Afcz-'= = 

fc=fco 



where Aj. G gl{m,¥), ko < k < ki, and k^ and ki are the smallest and the 
largest indices such that 7^ 0, respectively. We call the expression of A{z) 
in (12. ip the Laurent series of A{z). The number g = ki — ko + 1 is called the 
genus of A{z). If A{z) has genus g, then we can identify A{z) with the block 
matrix 

(2.2) A:={Ak,,...,Ao,...,Ak,), 
where 

(2.3) Ak = (ai^km+j), < i < m - 1,0 < j < m - l,ko < k < ki, 

and say that A{z) has size m x mg. We define the adjoint A{z) of A{z) by 

fcl — fco 

(2.4) Aiz) = A*iz-') = J2 Alz' = E 

fc=fco fc=— fcl 

where denotes the conjugate transpose of the complex matrix Ak. 

We are now ready to recall the definition of a biorthogonal wavelet matrix 
pair and the definition of a pseudoidentity matrix pair from 

Definition 2.1. A pair of m x mg matrices {L = {lij),R = (rij)) (as in 
) is said to be a biorthogonal wavelet matrix pair of rank m and genus g 
if the following conditions are satisfied: 

(2.5) L{z)R{z) = ml, 



mi 



and 



i j 



m, if i = 0, 

0, z/ 1 < 2 < m - 1, 



where Im is the m x m identity matrix. 

In the literature, the matrix L{z) is called the analysis matrix and the matrix 
R{z) is called the synthesis matrix of the biorthogonal wavelet pair. Condition 
(12.51) is called the quadratic or perfect reconstruction condition and condition 
(12.61) is called the linear condition. Note that (12. 6p , which is also referred to as 
the zeroth-order vanishing moment condition, is a necessary condition for the 
existence of scaling functions and wavelet functions. This is one of the main 
differences between wavelets and perfect reconstruction filter banks. 
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Definition 2.2. A matrix pair {C{z), D{z)) is called a pseudoidentity matrix 
pair if 

kc 

C{z) = J2 CkZ-\ D{z) = J2 DkZ-\ (fed < < K) 

(2.7) fc=0 k=ka 
C{Z)D{Z) = Irn, C{1) = D{1) = Im. 

We remark that it follows from the definition that 

(2.8) det{C{z)) = det{D{z)) = 1. 

For some other properties of the pseudo identity matrix pairs, please see 
Recall that a primitive paraunitary matrix is an m x 2m matrix of the form 

(2.9) V{z) = Im-vv* + vv*z-\ 

where v G F"^ is a unit column vector. Note that a primitive paraunitary 
matrix satisfies V{z)V{z) = 1^ 

The following theorem, due to Resnikoff, Tian, and Wells, reduces the fac- 
torization of a biorthogonal wavelet matrix pair to that of a pseudoidentity 
matrix pair. 

Theorem 2.1. (see [11, Thm. 3.4]) A pair of m x mg matrices {L,R) is a 
biorthogonal wavelet pair of rank m if and only if there exist primitive parau- 
nitary matrices Vi, . . . ,Vd, d> 0, such that 

L{z) = z~'^^V,{z) ■ ■ ■ V,{z)C{z) J ° ^ H, 

R{z) = z-''^V,{z) ■ ■ ■ Va{z)D{z) \ ^ H, 

where ko E Z, d = b — mko, b is the exponent o/det L{z), G G GL(m — 1, F), H 
is the canonical Haar matrix of rank m, and {Ci^z)., D{z)) is a pseudoidentity 
matrix pair such that the genus of C{z) is < g. 

We consider an example of pseudoidentity matrix pairs. 
Example 2.1. For any sequence of nonzero numbers 
(2.11) ai,a2,. . . ,afc e F, 

let < mi < ■ ■ ■ < rrik be k positive integers, let 



(2.10) 



(2.12) ao = ^a„ u{z) = J2 



k 

aiz 

i=l i=l 
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let 



(2.13) C7(z)=M"^° \^yfa, -a,\ 



and /ei 



1 + ao ao \ , / — ctj — ctj 



(2.14) - « i!l +E 



-ao 1 - ao / ^ \ ai a^ 

where cii is the complex conjugate ofai. Then 

C(l) = D{1) = h, C{z)D{z) = h, 

(2.15) 

det{C{z)) = det{D{z)) = 1. 

The first equality follows from the definition, and the second and the third 
equalities can be easily seen since 

1 — ao + u{z) ao — u{z) 
—ao + u{z) 1 + ao — u{z) 



(2.16) C{z 

has inverse 



(2.17) (C(^))-i= ( ^ + -ao + u{z) 

^ ^ V V ^ ao — ^(2;) 1 — ao + 'u(-2) 

and since D{z) = {C{z))~^. 
For instance, if we take A; = 2, ai = a2 = 1, and mi = l,m2 = 2, then 



(2.18) c(.)=( + ( i ill^'M J -1 1^" 

and 



(2.19) Diz) = { _\ ) + ( ) + ( ) z\ 



Remark 2.1. Since 

p-^o) (:: :::)'=o. 

these pseudoidentity matrix pairs provide counterexamples to the following con- 
jecture (Conjecture 1) in \1 Ij .• 

Conjecture: For m = 2 and C{z) and D{z) as in (2^, if C{z)D{z) = I2 
and p is the smallest positive integer such that C^^^p 7^ 0, then both Ck^-p and 
Dk^^p must be invertible matrices. 
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Remark 2.2. A weaker condition than the condition stated in the conjecture 
that guarantees the factorization of a pseudoidentity matrix pair was also given 
in [H] (see the remark after the proof of Thm. 4-4 [Hjy'- note that our 
examples do not satisfy this weaker condition either, since for the matrices in 
i\2.18\) the condition 

becomes 

But this last equality cannot be true, since if N is such that = 0, then 
h- N is mvertible with {h - N)'^ = 1^ + N. 

Remark 2.3. We note that the matrix in Ii2.18\) can be factored as 



-1 



(2.21) 

' ' 1 \ f 1 -1 \ / 1 -1 



I J V 1 -1 / V 1 -1 

Note also that the degrees do not add up. 



z-^ 



3. Complete factorization 

We now prove the following theorem. 

Theorem 3.1. Let {C{z), D{z)) be a pair of matrices of rank m > 2 defined 
by 

kc 

(3.1) Ciz) = J2CkZ-\ D{z) = Y,DkZ~^ h,k,>0. 

k=0 k=kj 

Then {C{z),D{z)) is a pseudoidentity matrix pair if and only if there exist a 
positive integer r and nilpotent matrices 

Ni G gl{m, ¥), = 0, l<i<r, 

such that 

C{z) = LnXz) ■ ■ ■Ln2{z)LnAz), D{z) = RnXz) ■ ■ ■ Rn2{z)Rni{z), 
with 

LnXz) =Im-N, + N,z-^\ RnXz) = Im + N*- Nlz''^ l<t<r. 
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Proof. It is clear that the conditions are sufficient. To prove the necessary part, 
note that since D{z) is just the adjoint of the inverse of C{z), we need only 
to prove the statement for C{z). To avoid writing the negative exponents, 
we make a change of variable (not essential) by letting t = z~^. We abuse 
notation and write 

kc 

(3.2) C(t) = J]Cfct^ 

Since the entries of C{t) are elements of the Euclidean domain F[t], by the 
Division algorithm, we can use elementary row and column operations of the 
form (type I elementary operations): 

add a multiple of a row or column to another row or column, 

to reduce C{t) to a diagonal matrix 



(3.3) C'{t) 



do 



\ dm / 

where 

(3.4) ^ rfi G F[t], l<i<m. 

This is possible since: (a) we can use elementary row and column operations 
to reduce C{t) to the given form (see for example [1, pp. 459-460]-here, 
the statement is actually weaker, since we do not require that di/di-i); (b) 
the operations that interchange two rows or columns can be obtained from 
type I elementary operations together with scalar multiplications by nonzero 
numbers to rows and columns; and, (c) in our diagonal form, we allow the 
elements di,l < i < m, to be defined up to nonzero scalar multiples (e.g. we 
do not care about negative signs), and thus multiplication by nonzero numbers 
to rows and columns are not needed here. 

Writing our reduction process as matrix multiplications, we have 

(3.5) C := C'{t) = P^it) ■ ■ ■ P,{t)C{t)P:{t) ■ ■ ■ P[{t), 

where the Pk{t), 1 < k < u, and Pl{t), 1 < / < s, are elementary matrices of 
the form: 

(3.6) Im + f{t)E,,, i^j, f{t)e¥[t], 

where Eij is the m x m matrix with 1 at the ijth entry and elsewhere. 

Since det(C(t)) = 1 (see (12.81) ). by taking determinants on both sides of 
(13.51) . we conclude that all (ij G F, 1 < i < m, that is, C is a constant diagonal 
matrix. 



8 FLORIAN M. SEBERT AND YI MING ZOU 

Let 

7^ pU)-f iPkit)r\ if l<k<U, 

[6J) ^fcW - I ^/(p/_j^))-l(c")-l, iiu + l<k<U + S. 

Then all Ek{t) are elementary matrices of the form (13.61) . and 

(3.8) C{t) = E,{t)---E^+s{t)C'. 
Consider elementary matrices of the form (13.61) . Let 

(3.9) Eit) = Im + fit)E,„ i^j, 
with 

(3.10) /(t) = a„r + a„_ir-^ + ■ ■ ■ + ao G ¥[t]. 

We can assume that f(t) ^ 0, since /(t) = gives the identity matrix, and 
factor E{t) as 

(3.11) E{t) = (/„ + (a„r - an)Eij) ■■■(/„ + {ait - ai)Eij){I^ + foo^ij), 
where a term exists only if the corresponding coefficient is not 0, and 

n 

(3.12) 6o = Xl«- 

Note that the factors in (13. lip involving t can be written as 

(3.13) Im - akEij + akEijt'' = Im-N + Nt\ 

where := a^Eij is nilpotent and A^^ = 0. If we conjugate the right hand 
side in (l37[3D by G G GL{m,¥), we have 

(3.14) G{In, -N + Nt^)G'^ =I^-N' + N't'', 

with (A^')^ = (GNG-^)'^ = 0. Therefore, upon plugging (KW) and (KW) into 
( 13.81) and moving all the constant matrices (e.g. the last term in ( 13. lip ) to the 
end, we obtain 

(3.15) G{t) = Lj^^it) ■ --LnMLnMQC' for some r > 0, 

where the factors L tv- {t) defined in (13131) or (13141) . and Q G GL{m,¥). 

Setting t = 1 in (I3.15P and noticing that Lj~^-{1) = /„, 1 < i < r, we have 
QG' = Im- Substituting back = t, we have proved the theorem. □ 

Remark 3.1. Since the factors are formed by nilpotent matrices, we should 
not expect that the degrees of the factors add up to the degree ofG{z) as \2.21\) 
shows. 

Remark 3.2. In general, the factorization depends on the reduction process 
and it is not unique. 
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Example 3.1. Here we consider a 3 x 3 example. Let 
C{z) = I 2^-2 + 2^-1-4 2^-2 H 



El = 



Eo — 



E. 



z 


-2 _ 


f - 2 ^-2 - 


f z 


1 







1 





1 















1 






1 ^ 








1 


1, 




z-2_^-i + 2 1 




1 











1 


2 








1 



Then 



E^E2EiC i^z^E^ — I^. 
Thus after moving the E^ 's to the right hand side and simplifying, we obtain 
C{z) = (74 - A^i + Niz-^){h -Ni + Niz-^){h -N2 + N2Z-^){h -N2 + N^z'''), 
where 

/ -1 -1 \ / -1 2 \ 

Ni^\ 2 2 , 7V2 = 2 -4 . 

\ 1 1 y \ 1 -2 y 

As a consequence of Thm. 3.1, we have a refinement of Thm. 2.1: 

Theorem 3.2. (biorthogonal wavelet factorization theorem) Form >2, a pair 
ofmx mg matrices {L, R) is a biorthogonal wavelet pair of rank m if and only 
if there exist primitive paraunitary matrices Vi, . . . , Vd, d > and nilpotent 

matrices 

Ni e gl{m, ¥), Nf ^0, l<i<r, 
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such that 



(3.16) 




H, 



where ko & 'Z, d — b — mko, b is the exponent o/det L{z), G e GL{m — 1, F), 
H is the canonical Haar matrix of rank m, and Lj^^.^Rj^., 1 < i < r, are as 
defined in Thm. 3.1 such that the genus of 



For a rank 2 pseudoidentity matrix pair (C(z), D{z)), the factorization can 
be described in a concrete way. We will see that the computational com- 
plexity in the factorization is basically determined by the operations involved 
in carrying out the Euclidean algorithm to determine the GCD of two poly- 
nomials with coefficients in F. Again, to avoid writing the negative expo- 
nents, we use the variable t instead of z in our discussion as before. Let 
a = a{t), b = b{t), c = c{t), d = d{t) e F[t] be such that 



Since C(l) = /2, ad ^ 0, and 6(1) = c(l) = 0. If one of b and c is 0, then 
det(C(t)) = 1 implies that a and d are constants and hence a — d—1. If, say 
c = 0, then we can write 



C{z) = LnX^) ■ ■■Lni{z) 



is < g. 



4. The rank 2 case 



(4.1) 




&(i) = (6X-M + --- + (M-&i), 



and factor C{i) as 



C{t)^LMSt)---LN,it)LM,{t). 



where 




1 < i < n. 
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So assume be 7^ 0. To fix our case of discussion, let deg(a) < deg(c) and 
apply the division algorithm repeatedly to get 

c = gia + ri, deg(ri) < deg(a), 
a = q2ri + r2, deg(r2) < deg(ri), 

(4.2) : 

rn-2 = qnTn-i + Tn, deg(r„) < deg(r„_i), 

Thus, by multiplying 

(4.3) {h - qn+iEuv) ■■■{h- q2Ei2){h - qxE^i), 

to Cif) from the left (where the index = 12 or 21 depends on the case), we 
obtain 

( X ) - ( I X 

Since the determinant of (14.31) is 1, we must have r„, d\ (or 61) G F. Multiplying 
(4.5) J2 - (rfi)-i6iEi2 or J2-(6i)-%E2i 

from the left again, we get a constant matrix: 

<«) ( 0^ : ) - ( : 0- 

Thus, letting the matrix that was actually used in the process from (14. 5 p be 

I2 - g{t)Ers, 

and the corresponding matrix from (14.61) be C we get 

(4.7) C{t) = [h + gi^2i)(/2 + ^2^12) ■■■{I2 + qn+iEuv){l2 + g{t)E,,)C'. 

Factoring each of the factors from the right hand side further as before if 
needed (see the proof of Thm. 3.1), we get the factorization desired. 

To summarize our discussion, we introduce some notation. For an integer 
m > 2, a 2 X 2 matrix Eij^i 7^ j, and a polynomial 

(4.8) /(t) = aX + --- + ait + ao GF[t], 
we define the factorization of Im + f{t)Eij as follows. Let 

(4.9) Lij-a^ = Irn " dkEij + t'^GkEij, 1 < k < Tl, Eijj(i) = Im + f{^)Eij, 

and set 
Then 

(4.11) Im + f{m,=L,,jit)E,,j^,y 
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Theorem 4.1. Let C(t) be defined as in ( f^.jp and assume ^4-'^ ^'^^ ( [^. 7| j 
hold. Then C{z) can he factored completely into a product of matrices of the 
form: 

h-N + z-'^N, = 0, 
by substituting for t in 

T T' T' T' 

where quA ^ k < n + 1, are as in (^T^; g = g(t),uv, rs, are as in (^T^; and 

-^21,(72 = -£^12,gi(l)-^21,g2-^12 

-^'l2,(?3 = -£'12,gi(l)-E'21,52(l)-^21,g2-^21,-g2(l)-^12 

are the conjugates of the corresponding Lij j. 

Proof. See the last paragraph of the proof of Thm. 3.1. □ 

5. Conclusion and Discussion 

The work of Resnikoff, Tian, and Wells [11] shows that any biorthogonal 
wavelet matrix pair can be decomposed into four components: an orthogonal 
component, a pseudoidentity matrix pair, an invertible matrix, and a con- 
stant matrix. Their work reduced the parametrization and factorization of 
a biorthogonal wavelet matrix to that of a pseudoidentity matrix pair. The 
main contribution of the current paper is the factorization theorem for a pseu- 
doidentity matrix pair of arbitrary rank proved in section 3. Our result implies 
in particularly in the case m = 2, that a slightly modified version of the factor- 
ization theorem stated in [11] is true in spite of the invalidness of the reduction 
procedure proposed there. We also provided several examples to explain our 
results. As the examples show, since the factorizations involve nilpotent ma- 
trices, in contrast to the orthogonal case [10], there is no uniqueness result 
even for the rank 2 case. Furthermore, in a factorization of a pseudoidentity 
pair, the degrees of the factors do not add up to the genus of the pair in gen- 
eral for the same reason. Hence the parametrization of these wavelets is not 
immediately clear from the factorization theorem. 

To explain this point in more detail, recall that [121 US] in the orthogo- 
nal case, the factorization of paraunitary matrices into products of primitive 
paraunitary matrices provides an onto map from 

^2m-l ^ ^2m-l ^ . . . x S^"^-^ xU(m) 

' V ' 

g factors 
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to the set of all paraunitary matrices of genus g, where U{m) is the set of 
unitary matrices (constant) and 

g2m^l = g W^lv^ = 1}. 

But in the biorthogonal case, the factorization does not lead to a similar result. 
More precisely, if (C(z), D{z)) is a pseudoidentity matrix pair of genus g, and 

C{z) = LnX^) ■ ■ ■ Ln2{z)Lnj{z), D{z) = RnXz) ■ ■ ■ Rn2{z)Rni{z), 

where the Ni, 1 < i < r, are nilpotent matrices such that A^^^ = 0, and 

LnXz) = In.-N, + N,z-'% RnAz) = Im + N*- N*z'% 1 < z < r, 

then we only have 

fci + ^2 H \- kr> g, 

and for some factorizations, strict inequality can hold. Thus, the parametriza- 
tion problem in the biorthogonal case is more complicated and needs further 
study, and the related geometry problems could offer solutions. These topics 
are currently under investigation. 

Finally, we would like to point out that though the factorization of a pseu- 
doidentity matrix pair provided in Thm. 3.1 is useful in analyzing the struc- 
ture of the biorthogonal wavelets such as parametrization, in applications, one 
may want to use the factorization provided by ( 13. 8 p and f l3.1ip directly. We 
state this alternative factorization as a theorem below. Note that because of 
the mix of the constant and nonconstant matrices in the factorization, the 
statement is not as neat theoretically. Let 6ij be 1 or according to whether 
i = j or not. 

Theorem 5.1. (lattice structure) The matrix C{z) of size m x gm is a pseudo 
identity matrix if and only if 

(5.1) Ciz) = Eriz)---E,iz)C', 

where the Es{z),l < s < r, are elementary matrices of the form 

(5.2) Ira + ihsfi - ^)asEij + asEijZ~^% i^ j, a^^O, < ks < g, 
and C is a constant diagonal matrix such that det(C") = 1 and 

( \ 



(5.3) 



n ^^(^) 

l<s<r , 



C — Im- 
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